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Abstract 

We propose the simplest model of teleparallel dark energy with purely a non-minimal coupling to 
gravity but no self-potential, a single model possessing various interesting features: simplicity, self- 
potential-free, the guaranteed late-time cosmic acceleration driven by the non-minimal coupling to 
gravity, tracker behavior of the dark energy equation of state at earlier times, a crossing of the 
phantom divide at a late time, and the existence of a finite-time future singularity. We find the 
analytic solutions of the dark-energy scalar field respectively in the radiation, matter, and dark 
energy dominated eras, thereby revealing the above features. We further illustrate possible cosmic 
evolution patterns and present the observational constraint of this model obtained by numerical 
analysis and data fitting. 
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I. INTRODUCTION 



The accelerating expansion of thepresent universe discovered in 1998 is one of the most 
important puzzles yet to be solved jl|. Its solution hopefully will lead us to a new revolu- 
tion in physics of this century. A positive cosmological constant, as a geometrical and/or 
an energy component of the universe, gives the simplest explanation that fits the current 
observational results. In addition, an energy source of anti-gravity dubbed "dark energy" 
and the modification of gravity provide two intriguing approaches to the solution 

A simple realization of dynamical dark energy is given by a scalar field minimally coupled 
to gravity, where the cosmic acceleration is driven by potential energy 3|, |4| or noncanonical 
kinetic energy j^J of the scalar field. For an alternative explanation from modified gravity 
one may introduce new degrees of freedom, modify the gravity action, or even change the 
formalism. A simple extension of general relativity (GR) in the first scenario is to introduce 
a scalar field non-minimally coupled to gravity, e.g., the scalar-tensor theory 6|. Models 
with non-minimal derivative couplings have also been proposed For the second scenario 
a straightforward modification is to invoke a nonlinear function of the Ricci scalar in the 
gravity action to replace the linear function in GR, i.e., the f(R) theory {§} that can be 
regarded as one of the scalar-tensor theories. As to the change of the formalism, an example 
is teleparallel gravity that is formulated with torsion and a curvatureless connection. 

Teleparallel gravity can be equivalent to GR when the gravity action invokes simply a 
linear function of the torsion scalar. To make an extension of teleparallel gravity as modified 
gravity, i.e. inequivalent to GR, one may follow the two scenarios raised above. A simple 
extension can be made by introducing a scalar field non-minimally coupled to teleparallel 
gravity. This can be regarded as a "scalar-teleparallel theory" of gravity, a modification of 
teleparallel gravity analogous to the scalar-tensor theory as a modification of GR. It has 



recently been proposed as an alternative dark energy model 9K12J and dubbed "teleparallel 
dark energy." For the other scenario a straightforward modification of the gravity action 
can be made by invoking; a nonlinear function of the torsion scalar, i.e., the f(T) theory (for 



a review, see 



m 



In this paper we focus on teleparallel dark energy. The distinct behavior of dark energy 
in this model has been investigated in {9), lo| where several different potential terms for the 
self-interaction of the scalar field are chosen. For lack of a guiding principle of choosing 



the self-potential, in this paper we propose the simplest model with no potential. That 
is, we investigate the "minimal" model of teleparallel dark energy where the scalar field 
is canonical, massless and noninteracting but non-minimally coupled to teleparallel gravity. 
With this simple model we attempt to manifest the effect of the non-minimal coupling on the 
cosmic acceleration, in contrast to the conventional models that attribute the acceleration 
to potential energy or noncanonical kinetic energy. 

In this model we find the analytic solutions of the scalar field in the radiation, matter, 
and dark energy dominated eras, respectively. These solutions largely facilitate the theo- 
retical studies on the dark energy behavior and help to reveal various interesting features 
of this model. In addition to the analytic solutions, we also numerically analyze this model 
and illustrate possible evolution patterns. As to be shown, the numerical solutions truly 
manifest the features read from the analytic solutions. We then present the observational 
constraint of this model obtained by the data fitting with the type la supernova (SNIa) 14 [ 



baryon acoustic oscillation (BAO) 15|, |l6[ , and cosmic microwave background (CMB) 17l.ll8| 
observational results. 



II. TELEPARALLEL DARK ENERGY WITHOUT SELF-POTENTIAL 

Teleparallel dark energy is played by a canonical scalar field non-minimally coupled to 
gravity in the framework of teleparallel gravity, which is formulated with the veirbein field 
eA(x^), the metric deduced from veirbein g^x) = r/AB^e^ , the curvatureless Weitzenbock 

w 

connection r£ = e^c^e^, and the corresponding torsion tensor 

w w 

1%, = - = e\ (d^ - d v e*) . (1) 

The action of teleparallel dark energy reads 

" T 1 

2 



S = I d A x e 



2 + ^(d^d^ + ^Tf) -V{<j>) + A 



(2) 



where e = det(e^) = yf—g, T is the torsion scalar defined by 

\rpppyrp , 1 



cj) is the scalar field, V{4>) the self-potential, C m the matter Lagrangian, and £ the non- 
minimal coupling constant. In this paper we investigate the simplest model free of self- 
potential, i.e., hereafter V{<p) — 0- 



For a flat, homogeneous and isotropic background where the space-time is described by 
the flat Robertson- Walker metric ds 2 = dt 2 — a 2 (t)5ijdx l dx : > and the scalar field is simply 
time-dependent, the scalar field and the gravitational field equations |tj read 



where the matter energy density p m oc a 3 , the radiation energy density p r oc a 4 , and the 
energy density and pressure of the scalar field are given by 



We note that for £ < 0, the non-minimal coupling may result in negative pressure, thereby 
anti-gravity, and meanwhile a positive energy density In particular, the contribution from 
the term 3£H 2 <j) 2 entails the cosmo logical-constant equation of state w = —1 and the other 
term 2£d(H<j) 2 ) / dt provides no energy density but negative pressure when Hip 2 increases with 
time. Henceforth we consider the case where £ < in order to guarantee the positiveness of 
the energy density. We set the present scale factor a = 1 for convenience and without loss 
of generality. 

III. ANALYTIC SOLUTIONS 

Here we present the analytic solutions of the scalar field in two cases: (a) the case where 
H oc 1/t, including the radiation-dominated (RD) and matter-dominated (MD) eras, and 
(b) the scalar-field-dominated (SD) era. 

(a) H = a/t, i.e. a(t) oc t a , with constant a: 



4> + 3H<p + 6£H 2 <p = , 




(4) 



P0 = ^ 2 "3£# 2 2 , 

p^ = \ft + ^H 2 <p 2 + 2ij t (H<p 2 ). 



(5) 



<j>{t) = dt h +C 2 t l \ 



(6) 



where are constants and 



±^(3^ - l) 2 - 24£a 2 - (3a - 1) 



(7) 
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In RD, a = 1/2; in MD, a = 2/3. For £ < 0, the power-index l\ is positive and I2 negative, 
corresponding to an increasing mode and a decreasing mode, respectively. 

Hereafter we take into consideration only the increasing mode, i.e., <p(t) = Cit l1 . The 
energy density, pressure, and equation of state corresponding to this solution are given by 



[1 2 J2 + 3£a 2 + 2£a(2h - 1)] Cft 2 ^ , 



V4> 



4£a 



2(l-/i) 
3a 



Consequently, one obtains 

RD: wj, = - (2 
MD: = - (l 



24£ 



oc a 



-5+V1-245 



32f/3 



oc a 



(-9+V9-96£)/2 



(8) 

(9) 
(10) 



To sketch the behavior of teleparallel dark energy in RD and MD, in Tabled] we consider 
four different £ as the examples for demonstration. The £ values under consideration are 
within the range [—1, 0) because the case where £ < — 1 is rather disfavored by observational 
results about the cosmic expansion. The four values of £ are so chosen that in RD or MD the 
dark energy density behaves like the familiar energy sources such as radiation, matter, or a 
cosmological constant: (i) When £ = 0~, teleparallel dark energy tracks the dominant energy 
components respectively in RD and MD, a behavior similar to the tracker quintessence 19] 

(ii) When £ = —1/8, in RD p^ behaves like matter and in MD it drops slower than matter. 

(iii) When £ = —3/4, in RD p^ decreases rather slowly and in MD it is constant, (iv) When 
£ = — 1, in RD p^ is constant and in MD it increases as phantom. 

TABLE I. Explicit solutions of the scalar field in RD and MD with different values of £. 





(i) 0- 


(ii) -1/8 


(iii) -3/4 


(iv) -1 


Era 


RD 


MD 


RD 


MD 


RD 


MD 


RD 


MD 


h 


-3£ 


-8£/3 


+1/4 


+0.264 


0.840 


+1 


+1 


+1.208 




1/3 + 4£ 


-8£/3 





-0.264 


-0.786 


-1 


-1 


-1.208 


P<t> oc 


a" 4 a- 12 « 




a" 3 


a" 2 - 21 


a -0.641 


Constant 


Constant 


a 0.623 



The analytic solution, as well as the examples in Table HJ shows several interesting fea- 
tures: The dark energy equation of state is a constant simply determined by the sole 



model parameter £ but insensitive to the initial condition of the scalar field, i.e., a tracker 
behavior of the dark energy equation of stated For £ < 0, w<j> is smaller than that of the 
dominant energy component, w a = — 1 + 2/(3a), at all times in RD and MD. Accordingly, 
the dark energy density p^ decreases slower than the dominant energy density p a and there- 
fore the late-time domination of dark energy (SD) is guaranteed. For larger |£|, is smaller 
and p<t>/ p a increases faster, and therefore SD comes sooner for a fixed initial condition. 

The analytic solutions in RD and MD describe the universe at earlier times. At later 
times, the universe is dominated by both matter and dark energy (MSD), e.g., in the present 
epoch, and then dominated solely by dark energy (SD) in the future. In the following we 
attempt to analytically integrate the field equations in both of the later eras (MSD and SD). 

(b) Scalar field domination: 

For the universe dominated by matter and dark energy, we deduce an integrable equation 
from Eq. (TJ): 

| [F{^H] = ^p m a 3 = \tf$U , 

F{<f>) = 1 + K 2 i4>\ (11) 

where Hq and Q m o are the Hubble constant and the matter density fraction at the present 
time, respectively. Its solution gives a relation: 

F(<f))a 3 H = ^H 2 Q m0 t + C 3 , (12) 

where C 3 is an integration constant. This can hardly be further integrated analytically. 

For the universe dominated by dark energy solely, Eq. ([121 gives a simpler relation with 
^mo = 0: 

F((f))a 3 H = C 3 . (13) 
From this relation and the field equations we obtain the second integrable equation, 

J-(J*X = ± (14) 

from which one can straightforwardly obtain the solution, 

/ / v ,sin6 
<p(a) = ±- 



9{a) = ^6£lna + C 4 , (15) 



1 This is different from tracker quintessence [l9( where it is the scalar field (thereby the dark energy density 
and its equation of state) that possesses the tracker behavior. 
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and then the equation of state, 



Wf = -l- >/ = 32l73 tan 9, (16) 

where C4 is an integration constant. 

This solution predicts a future singularity where both the expansion rate H and the dark 
energy equation of state w$ go to infinity It is one kind of finite-time future singularities, 
i.e. "type III" singularity classified in Ref. [20 1 . The singularity occurs when 9 — (n + 



for integer n, i.e., = ±l/y — k 2 ^ and F((f)) = 0. Between two singularities, e.g., for 
— 7r/2 < 9 < 7r/2, <fi is monotonic and monotonically decreases from +00 to —00 along 
with the expansion, experiencing the crossing of the phantom divide = — 1 at 9 = 0. 

IV. POSSIBLE EVOLUTION PATTERNS AND DATA FITTING 

Here we illustrate the evolution patterns of the dark energy equation of state for 
a wide range of initial conditions by numerically solving the field equations 01]). These 
numerical solutions will manifest the interesting features we have read from the analytic 
solutions, particularly the feature that approaches the constant tracker values in RD 
and MD, leaves the tracker and starts decreasing when the universe is leaving MD for SD, 
keeps decreasing and crosses the phantom divide in SD, and goes to negative infinity at a 
finite time eventually. In addition, we confront this self-potential-free teleparallel dark energy 
model with the observational data that are relevant to the cosmic expansion, including SNIa 



BAO 



15 



16| and CMB 



17 



thereby obtaining the constraint on the sole model 
parameter, i.e. the non-minimal coupling constant £, as well as the cosmological parameters. 

For illustrating the evolution patterns of numerically, we consider the case where 
£ = —0.35 and show io^(log(a)) in Fig. [T]for different initial conditions. For the same case we 
also show the present values of (Q m , w^) in Fig. |2j We consider two sets of initial conditions 
at the initial time log(a) ~ —8.69 (i.e. ln(a) = —20): (i) n<f) = and K(j)' = d(n(j))/d\na G 
[10" 50 ,2.5 x 10~ 8 ] (solid line), and (ii) nfl = and k0 e [10^ 50 , 1.2 x 10~ 8 ] (dashed line). 
For these two kinds of initial conditions, initially w,p is unity and 1/3, respectively. Then, as 
shown in Fig.fU it quickly evolves to the RD tracker value, —0.36, along with the domination 
of the increasing mode of 0(t)-field over the decreasing mode. It reaches the tracker value 
around log(a) = —8. When the universe gradually leaves RD and enters MD, gradually 



decreases to the MD tracker value, —0.58. Later, when teleparallel dark energy starts to 
dominate, leaves the MD tracker and decreases monotonically. Finally, crosses the 
phantom divide (w$ = — 1) and goes to the singularity as expected. 
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FIG. 1. Evolution of the dark energy equation of state, w ( p(log(a)), with £ = —0.35. The black 
area in SD is formed by the trajectories of w ( j ) (log(a)) w.r.t. a wide range of initial conditions. 
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FIG. 2. The present values of (£l m ,W(fi) for a wide range of initial conditions with £ = —0.35. The 
curve shows a relation between these two quantities, as a consequence of the tracker behavior of 
for a flat universe with given £. 

Because of the tracker behavior of w^, the evolution of dark energy is simply determined 
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by two parameters, the model parameter £ and the integration constant C\ appearing in 
Eq. ([6]), where C\ may be regarded as the only important initial-condition parameter. Ac- 
cordingly, for a flat universe with given £, the present value of and the present matter 
density fraction Q m0 are related by a single parameter C\. This relation is shown in Fig. |2J 
where both Q m0 and w^q change with the initial condition, i.e. the value of C\. The present 
value of has an upper bound that depends on £ and is set by the tracker value of in 
MD. Note that the two different sets of initial conditions give the same relation curve. 



Following the same procedure of data fitting in Ref. 10J, we obtain the observational 
constraint of our model. Figure [3] shows the la-3cr confidence regions in (f2 m o, £) (left panel) 
and (n, m o, w^o) (right panel) obtained from the SNIa (blue), BAO (green), CMB (red), and 
the combined (black) data. The best fit to the combined data locates at £ ~ —0.35 and 
f2 m0 ~ 0.28. The three data sets do not give a good concordance region, indicating the 
imperfection in simultaneously fitting these data sets. This is possibly due to the feature of 
this model that p^ may not be negligible at early times and accordingly may significantly 
affect the CMB result. We also remark that in contrast the data can be fitted well in the 
teleparallel dark energy models with potentials in Refs. [si, [li]]. This is expected because the 
models with potentials are more flexible than our simple potential-free model when fitting 
data. 
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FIG. 3. The la-3a confidence regions in (fi m o>£) (left panel) and (n m o,u^o) (right panel) 
obtained from the SNIa (blue), BAO (green), CMB (red), and the combined (black) data. 
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V. SUMMARY 



We propose the simplest model of teleparallel dark energy where dark energy is played by a 
scalar field that is canonical, massless and noninteracting but purely non-minimally coupled 
to teleparallel gravity. We have found the analytic solutions of the scalar field, thereby 
obtaining the analytic behaviors of the dark energy equation of state, in the radiation, 
matter, and dark energy dominated eras, respectively. With the analytic solutions we have 
shown various interesting features of this simple model: 

• The cosmic acceleration is driven by the non-minimal coupling to gravity, but not 
by potential energy or noncanonical kinetic energy as invoked in the conventional 
scalar-field dark energy models. 

• The domination of dark energy, the crossing of the phantom divide, and therefore the 
occurrence of the cosmic acceleration at late times are destined. 

• In the radiation and matter dominated eras the dark energy equation of state is roughly 
a constant and has tracker behavior, i.e. being insensitive to the initial condition. 

• In the dark energy dominated era the dark energy equation of state decreases mono- 
tonically and crosses the phantom divide at a late time; its value strongly depends on 
the initial condition of the scalar field. 

• The universe will meet a future singularity with the expansion rate going to infinity 
within a finite time. 

• This model has simply one free parameter, the non-minimal coupling constant £, for 
a spatially flat FLRW universe with the present matter density Q m o fixed. 

We have fitted the model parameter £ with the observational data of SNIa, BAO and 
CMB in cosmology. The concordance region for all these three data sets is only at the 3a 
level. It may indicate a slight incapability for this model to describe the early universe. 
This requires further examination. For further investigations, in addition to cosmological 
observations, the local gravity tests are expected to give significant constraints on teleparallel 
dark energy and therefore should be worth detailed studies. 
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